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This paper aims to get an optimal high voltage magnetic flux leakage 
transformer design of a three-phase shell type. Optimal design of transformer 
requires determination of design variables to optimize a particular objective 
and satisfying a set of constraints. The objective function is to minimize 
the total mass and reduce the volume of the transformer. This function 
depends on inputs, which are divided into optimization variables. 
Each optimization variable varies within a certain interval thus defining 
a global search space. It is within this space that we seek the optimal 
solution. The constraints: maximum and average current of magnetron anode 
are part of the problem in order to limit the overall search space. The results 
obtained indicate that the method has provided a global optimum. 
The computation time and cost of active material are much reduced 
compared with the conventional design results. 
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1. INTRODUCTION 

Optimization presents the way to find out the best possible optimum design, and at the same time 
satisfies all limitations or restraints inposed on its performance. Several techniques have progressed through 
centuries to deal with the optimal solution. Some of them are traditionally based on mathematical formulation 
and others are non-traditional based on soft-conputing techniques. 

As far back as the beginning of the 20th century, early research in transformer design attenpted to 
reduce much of this judgment in favor of mathematical relationships [1, 2]. The manufacturers started to 
research optimization methods [3, 4]. They sought to reduce much of this judgment with analytical formulas 
using a conputer was a pioneer by [2, 5]. Later [6-9] suggested a method for getting an optimized design of 
distribution transformer. Several other techniques were also used [10-13]. 

Looking for the optimal solution requires many iterations. It is difficult to accelerate the same 
through sinple calculations using a calculator. To get an optimal solution option must be made to a digital 
conputer. It can perform calculations at an extremely high condition and it has got a large amount of 
memory. In fact, it is the advantage of the conputer, which has revolutionized the field of optimization. 

The magnetic flux leakage transformer (MFLT) is the most inportant and costly components used 
in power supply systems. A good transformer design satisfies certain functions and requirements, such as 
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powering many magnetrons by keeping the same transformer with respecting the constraint inposed without 
damaging them. We can satisfy these requirements with various designs [14, 15]. The aim of this study is to 
find the most economical choice and arrive at an optimum design respecting the limitations inposed 
by the constraint functions. As a result is to reduce the volume, ensure better utilization of material, it saves 
space, and reduce loss, consequently get the lowest cost unit. In the case of the MFLT, we considered several 
aspects ofthe transformer design such as core shape, size, properties, copperwire turn, etc. 

Minimization of the volume and the weight depend on inputs that are divided into optimization 
variables. These variables are the size of the magnetic circuit characterized by the width of outer limb (a), 
the number of secondary turns (n 2 ), the size of each gap between the two shunts and the magnetic circuit (e) 
and the number of stacked sheets (ns). 


2. GEOMETRY OPTIMIZATION OF THE MAGNETIC FLEX LEAKAGE TRANSFORMER 

The cross section for shell type transformer is a rectangular shape as shown in Figure 1, 
which is a typical three-phase three-limbed core structure used for stable power supply applications. The core 
is made by stacking varnished laminations of silicon steel SF 19 . Either copper is used as conductor material. 
The primary and secondary windings are wound in the center of the core columns. Different from 
the ordinary transformer, the flux leakage flowing in the magnetic shunts could not be ignored [14, 16]. 
There are two vertical shunts for each phase conposed of silicon steel sheets in the center. The inportance of 
steel sheets is to provide a magnetic flux path, which laminated along the direction to reduce the iron loss. 
The air gaps, one at each end of the shunts, provide for different flux densities in the primary and secondary 
portions of the center leg. 

Figure 2 presents the equivalent model quadruple in tt of three-limb three-phase power supply 
treated and developed in previous papers [16]. Where uij and U 2 j are the voltages at the transformers 
terminals, m, i’lj, r’lj, L’pj and n 2 , Ej, r 2 j, Lsj are the number of turns, the currents through the coil, 
the resistances and the leakage inductances of the primary and the secondary windings, respectively. 
L’shj, i’sh, i’pj and hj are leakage inductances of the shunt, the currents through the shunt, primary and 
secondary inductance, respectively. Subscript] denotes the phase. 




Figure 1. Investigated shapes Figure 2. The equivalent electrical model of the three-phase 

of three-phase magnetic flux transformer power supply one magnetron per phase [17, 18] 

leakage transformer (Geometry 
and dimensions) 
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The obtained equivalent model nonlinear characteristics were inplemented in MATLAB Simulink 
to simulate the transformer response of current and voltage waveforms recorded during simulation 
tests [16-18]. The results confirm that the reference model is available [14, 15, 17, 18] and can use it in order 
to get more results that are satisfying. This is the objective in the following sections. The basic design 
variables in this optimization method are: X = {x^,X 2 ,X 3 ,X 4 }, where xi is the width of the uncoiled 
core a, is the number of turns at the secondary n 2 , xs is the size of each gap e and X 4 is the number of 
stacked plates ns. 


3. MFLT DESIGN OPTIMIZATION PROBLEM 

The aim of transformer design optimization is usually to optimize an objective function, 
which is subjected to several constraints. Among the various objective functions, the frequently used 
objective functions are the minimization of volume, total mass, main material cost, space, congestion and 
maximization of transformer nominal power. Under parameter constraints, the problem is e:?q)ressed 
as follows: 

Minimize (or maximize) /(X) 

X} < X KXj^, i=l...n 


Where X presents the vector of transformer parameter, f(x) is the objective function to be minimized or 
maximized, Xj andX^ are lower andupper bound of j-th parameter Xj. 


3.1. Optimizing techniques 

This subsection proposes a new optimization approach [7, 10, 19] for this special HV magnetic flux 
leakage transformer. The objective function is to minimize the total mass of copper and core material, 
by minimizing at the same time aU geometric parameters of the transformer. Searching for the optimum 
design using an appropriate mathematical technique. The process of optimization initiates by choosing 
the optimization variables and specifying the design constraints. The chosen variables are to be defined and 
their bounds determined. Then the objective function is to be formulated. 

The principle of optimization technique proposed is to perform a set of simulations. We use 
essentially the reference case [14] of the equivalent model to study the sensitivity of the magnetron current, 
due to variations of one or more geometric parameters. Eigure 3 shows the results of variables variation of 
the average and maximum current value at the nominal mode. We study the influence of each parameter by 
varying only one; the other parameters remained unchanged and identical at those specified in 
the reference case. It concludes that the variation of these parameters changes the electrical functioning of 
the high voltage power supply. In each simulation, we observe the waveforms of different electrical 
parameters of the circuit HV especially those giving the shape of the current magnetrons. 


3.2. The objective function formulation 

Each optimization design variables vary within a boundary condition. In order to frame the objective 
function in terms of the design variables, only solutions that verify the constraints defined. This problem can 
be defined by a triplet (X, D, C) where: 

X = {xi,X 2 ,X 3 ,X 4 } is the set of variables of the problem, where x^ represents the width of the unworn 
core (a), X 2 indicates the number of turns at the secondary (n 2 ), X 3 is the size of each gap (e) and X 4 
corresponds to the number of stacked sheets (ns) 

The fields of definition: D = [D^, D 2 , D^, D^} is the set of domains of the variables, for aU /c G [l; 4] we 
have Xj,eDj,, thus = [l5:2.5: 25], = [2050 :200: 2450], = [0.45:0.25:1.05] and 

= [10:2:18]. 

Each of these variables is assumed to be continuous. The following constraints C = { 61 , 62 , 63 } 
are inposed on the design problem 


/ = fi(.X^,X2,Xo,X^) 
Im=Imin=f2 0) 

— ^aveg ~ /s 

Vol = f^ix^,X 2 ,X^) 


- C, : < I^< 

r . jmin ^ j ^ jmax 

C 2 • ^ h - hs 

63 : 61,62 and min(VoZ) 
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With: 

The function is provided by Simulink of the magnetron current (A) for the equivalent global model. 
The function ^ calculates the minimum value of the magnetron current in the interval [0.6, 0.64 (s)]. 

The function calculates the average value of the magnetron current in the interval [0.6, 0.64 (s)]. 

Ci: Constraint on the minimum current (A) which guarantees the operation of the magnetron at full 
power and whose values should belong to the interval [/ I 1.2 -0.80]. 

C 2 : Constraint on the average current (mA) which guarantees the operation of the magnetron at fuU 
power and whose values should belong to the interval [/ I 300 -200]. 

C 3 : Conpound constraint which inposes the satisfaction of Ci and C 2 and which seeks the states or 
the state having the minimal volume. 






Figure 3. Influence of study design variables on the behavior of magnetron load maximum and 

average currents 


The function calculates the volume total (cm^) of the transformer conposed of the copper and 
tank core volume as shown in the (1). 

Vxot = Veu + Vpe (1) 

Transformer tank volume can be obtained as defined with the following formula 

Vpe = (A ■ C + 6 ■ h ■ (F - e) - 6 ■ E • F) ■ B (2) 
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The (2) can be rewritten in terms of design variables (xi, andx 4 ) as following 

Vpe = [72x2 + 3 . X 3 . (x^ - 2 X 4 )] ■ 60 (3) 

The volume of copper wire in a transformer depends upon the mean winding length Ln of primary and 
secondary windings, the total number of turns and cross -sectionalarea. It is given by 

Vcu = Lpj • Spj + Ljj • Sjj (4) 

Calculating the length of primary and secondary turns for each phase 

Lpj = ■ n^ and L^j = Lj^j ■ n 2 

Calculating the mean length per turn (Lng) for both primary and secondary coils 

= 2 ■ B + 2 ■ D + 4 ■ F (5) 

The volume total of copper wire in a transformer 

Vcu = Lmj ■ (ni ■ Spj+n 2 -Ssj) (6) 

The ( 6 ) could write also according to xi and 

Vcu = 3 ■ (120 + 8xi)(424.8 + 0.20x2) (7) 

EinaUy, from the (3) and. (7), we get the volume total (cm^) of the transformer: 

f 4 (xi,X 2 ,X 4 ) = [72Xi 2 + 3 • Xg • (Xi - 2 X 4 )] • 60 + 3 • (120 + 8Xi)(424.8 + O. 2 OX 2 ) ( 8 ) 

Calculating the total weight of copper in the transformer, mcu 

^cu “ Pcu ' ^cu (^) 

^cu = 3 ■ ■ (120 + 8xi)(424.8 + 0 . 20 x 2 ) (10) 

Calculating the weight of the iron core in the transformer 

Hlpe = pFe ' [’72Xi + 3 ■ X 3 ■ (x^ - 2 X 4 )] ■ 60 (11) 

The total weight is taken as the objective function, which is influenced by the design variables X and 
constraints.The function / is written like the following: 

mpe + = pFe^Fe + Pcu^Cu ( 12 ) 

The optimizing function has been conputed and the expression is given below 

/(X) = ■ [72x1 + 3 x 3 ■ (Xi - 2 x 4 )] ■ 60 + Pcu ■ [3 ■ (120 + 8Xi)(424.8 + O. 2 OX 2 )] (13) 

Where Pcu = 8940 is the mass density of copper in [Kg/w?] andpp^ = 7650 is the mass density of core 
material in [Kg/m^]. 

3.3. Procedure for optimization 

In order to get an optimal solution, the first step needs to formulate the problem, choose the design 
variables, define the constraints and precise the objective function. Minimum (or Maximum) bounds 
are inposed on the design variables. This subsection describes the method for the optimal design of 
a three-phase shell-type distribution MELT using the method of sequential programming (SQP) [20-25]. 
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The use of this method is successful in optimization of dimensioning and shape optimization. This algorithm 
is powerful and effective in the nonlinear programming, attenpts to resolve the program directly instead of 
transforming it into a sequence problem of minimization without constraints, which makes this algorithm 
differs conpared to other methods (method of optimization without constraints). The advantage of the SQP 
method is that it can be manipulated in MATLAB by using the function "fmincon" in the toolbox of 
MATLAB. This algorithm minimizes a given objective function respecting the constraints determined by 
the user, where the objective function, defined the total volume of the transformer with shunts, 
is in the following form The algorithm can be explained in the form of pseudo-code as follows: 

Algorithm: 

a. Read transformer data, independent variables, constrain and define of domains of the variables. 

b. Set the loop conditions for each variable for aU Xj^E Dj^. 

c. Check constraint (C 3 ) which inposes the constraints of (Ci) and (C 2 ) and which seeks the states or 
the state having the minimal volume. 

d. Accept or reject each point, check if the C 3 is verifier and go to next step. Otherwise, go to step (b). 

e. Formulate of the objective function and calculate the initial volume. 

Once this analysis is completed, the decision is making the ability of the conputer to get the best 
solution with respecting the constraints inposed in order to obtain the lower volume and mass. 

Performance characteristics obtained from the simulation test has been listed in Table 1. 


Table 1. Conparison of results obtained using the conventional design of three-phase 


transformer performances 


Solution 


Sdata 

Si 

S2 

S3 

S4 

Design 
variables X 

Xl 

25 

22.5 

20 

17.5 

15 

X2 

2400 

2400 

2750 

2050 

2400 

[an2n3 e] 

X3 

18 

18 

10 

18 

18 

X4 

0.75 

0.9 

0.9 

0.45 

0.75 


Phase 1 

1.0418 

0.932 

0.931 

0.829 

0.785 

Imax (A) 

Phase 2 

1.0469 

0.933 

0.932 

0.833 

0.785 


Phase 3 

1,0459 

0.931 

0.931 

0.829 

0.788 


Phase 1 

298,1 

265.1 

280.1 

201.3 

205.2 

Imean (mA) 

Phase 2 

297,9 

265.1 

280.0 

201.5 

205.4 


Phase 3 

298,3 

265.3 

280.3 

201.5 

205.3 


Phase 1 

1240,9 

1092.9 

1157.6 

820.4 

835.5 

Pmean (W) 

Phase 2 

1239,9 

1092.8 

1157.4 

821.3 

836.2 


Phase 3 

1241,8 

1093.6 

1158.4 

820.7 

835.8 

Volume 

(cm^) 

Vcu 

868.6 

814.3 

818.8 

651.1 

651.4 

Vpe 

2776.1 

2254 

1760.7 

1376.7 

1015.7 

Vxot 

3644.7 

3068.3 

2027.9 

2027.9 

1667.2 


mcu 

7.76 

7.26 

7.32 

5.82 

5.81 

Mass (Kg) 

mpe 

21.24 

17.24 

13.47 

10.53 

7.78 


f(x) 

29 

24.5 

20.79 

16.35 

16.59 


It is supposed that we have a calculation for a three-phase 1650VA transformer. The input voltage 
source is 50Hz, 220 V sinusoidal AC voltage with a +10% regulating range. The secondary voltage is 2230V. 
The transformer has 3 phases, 3 limbs; the silicon steel used for the iron core material is SF 19 . A voltage 
doubler capacitor C=0.9pF and high voltage rectifier diode Di. The main program, objective function 
calculation, and optimization algorithm were coded in MATLAB 2014a, including the inplementation of 
the equivalent electrical model of the three-phase transformer power supply. 

3.4. Dimensions of shell type transformer 

During this work, we have taken the manufacturers' data as a reference. The following geometrical 
dimensions of the shell type, both optimized and reference transformer have been shown in Figures 4 and 5. 
Core design includes the window height, window width, and the outer and central limbs. The winding design 
includes the area of cross section and number of turns. The geometry design and dimension details of 
the optimal transformer are given below. The transformer has been obtained by precisely considering 
he resulted dimensions calculated from the practice described in sections above. Measurements of 
he dimensions obtained from the calculations have been enlisted in Table 2. 
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A=375iiim 

B=60mm 

C=150iiim 

D=50mm 

E^75mm 

F=25mm 

hwp^w^33nim 

h=9mm 

e=0.75mm 

np=240 

ns=2400 

Vcu=868^^ 

VFe=2776.1cm^ 

Vtot=3644.7cmi’ 



A=300nim 

B=60mm 

C=120nim 

D=40mm 

E>:60mm 

F=20mm 

h wp=iiw^275iiim 

h=5mm 

e=0.9mm 

np=240 

ns=2750 

Vcu= 818 ^^ 

VFe=1760.7cm^ 

Vtot=2579.6cm^ 


Figure 4. Dimensions of shell type magnetic shunt 
flux leakage transformer reference case 


Figure 5. Dimensions of shell type magnetic flux 
leakage transformer optimized case 


Table 2. Conparison with a manufacturer’s data and optimized transformer 




Resulted Measurements 

Dimensions 


Manufacturers' 

data 

Optimal Design 


Width of outer limb (D/2) 

2.5 cm 

2 cm 

Core Design 

Depth of core (B) 

6 cm 

6 cm 

Width of central limb (D) 

5 cm 

4 cm 


Central limb cross section (2ab) 

30 cm^ 

24 cm^ 


Height of window (E) 

7.5 cm 

6 cm 

Design 

Width of window F) 

2.5 cm 

2 cm 

Area of window (Aw) 

18.75 cm^ 

12 cm^ 

Overall 

Overall Width (A) 

15 cm 

12 cm 

Frame Design 

Overall Height (C) 

37.5 cm 

30 cm 


Cross section primaiy (Sp) 

1.77 

1.77 


Cross section secondary (Ss) 

0.20 mm^ 

0.20 mm^ 


Diameter primary (cpp) 

1.5 mm 

1.5 mm 


Diameter secondaiy (cps) 

0.5 mm 

0.5 mm 

Windings 

Design 

Secondary Resistance (Rs) 

650 

64.68 0 

Primary Resistance (Rp) 

1000 

89.3 0 

Primary Wire Length (Lp) 

71.68m 

62.72 m 


Secondary Wire Length (Ls) 

768 m 

770 m 


Mean Wire Length (Lm) 

320 mm 

280 mm 


No. of primary turns (np) 

224 

224 


No. of secondary turns (ns) 

2400 

2750 

Shunt Design 

Height of shunt (h) 

9 mm 

5 mm 

No. of stacked sheets (ns) 

18 

10 


4. COMPUTATIONAL RESULTS AND DISCUSSION 

Simulations have been conducted in order to verify the proposed optimal transformer model. 
The work started with study the influence of various geometrical parameters of the transformer. Then we 
observe their influence on the magnetron currents. The proposed algorithm is selected to study. The high 
voltage transformer response of the solution chosen is listed during the simulation. Therefore, Statistical 
results and design variable values for the best solutions can be highlighted in Table 1. 

The algorithm takes into account many variations in design variables. These variations allow 
the investigation of a candidate solution. For each one of the candidate solutions, it is verified whether 
the constraints are satisfied, and if they are satisfied, the volume and mass is calculated and the solution 
is considered as acceptable. Finally, among the acceptable solutions, the solution S 2 with the minimum 
manufacturing cost is selected, which presents the optimum transformer. S 2 allows the best conpromise 
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between the total volume of the transformer and the operation of the magnetron. As a result, the objective 
function calculations would be reduced by 28.3%. Another 6.71% reduction is obtained for lost energy units. 
Based on these values of variables. The solution S 4 presents a minimum volume, but it does not allow 
the functioning of the magnetron in fuU power during its operating life. 

Using the decision variables of the solution S 2 , We simulated under MATLAB Simulink 
the electrical behavior of the HV optimum transformer for a magnetron [25]. The waveforms corresponding 
to this solution are shown in Figure 6. The waveforms obtained shows that these results are in perfect 
accordance with those obtained by the reference case, and respect the criteria recommended by 
the manufacturer for each phase (Imax<1.2 A and Inean < 300 mA.) Each magnetron into a phase that operates 
at rated speed (220V and 50Hz in the primary side), the electrical signals in the diode, the capacitor, 
the magnetron, and the secondary side have the same shape as those of a conventional single-phase power 
supply for one magnetron. These signals are periodic but not sinusoidal and they are phase shifted 120°. 


Magnetron 



Time (s> 


Secondary 



Time (s) 



0.6 O.&I D.62 0.54 


Tjme 



TInoe <s) 




6000 
4000 
. 2000 

-2000 

-4000 


0.61 0.62 0.63 0.64 

Time <S) 


Figure 6. Simulation curves voltages and currents of the optimized transformer conpared with 

the Manufacturers'data transformer 


It is noted that the maximum allowable value of the anplitude of each magnetron’s current in 
the optimized transformer does not exceed the acceptable limit (<1.2A), which conplies with the constraints 
inposed by the manufacturer. It ensures the correct functioning of the magnetron with a reasonable average 
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current of 300mA without exceeding the recommended peak current. Qven the above, the current 
stabilization process in each magnetron is completely ensured, which conpletely protects this 
microwave tube. 

The magnetron currents regulating the process in each phase of this optimal solution was verified. 
While observing the stability of the current’s variations in each magnetron with respect to the variations of 
the primary voltage of +10% around the nominal voltage of 220 Volts. Eigure 7 shows the waveforms of 
each magnetron’s current corresponding to the respective values of 200V and 240V on the primary voltage. 



Time (s) 


0.4 

0.2 


Time (s) 


Magnetron 


Magnetron 


Eigure 7. The voltage and current of magnetron anode for each phase with different voltage 

source values 


5. CONCLUSION 

In this paper, the detailed technique of a transformer optimization model of three-phase, shell-type 
transformers are presented, which can make a conpetitive solution. The attenpt of this paper is to optimize 
the volume of this transformer and establish that this technique is capable of giving more solution. Thus 
is a viable tool for obtaining the optimal design of the three-phase transformer. Based on the value of 
the design variables found out by the algorithm, the dimensions of new optimized transformer selected and 
the performance variables have been calculated. Eor the validation of the model, simulation results 
are conpared with published measurements. Although the paper is written for three-phase transformers, 
the developed model is suitable for the simulation of any other configuration of three-phase magnetic flux 
leakage transformer with several magnetrons by phase. 
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